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Abstract 

In this paper we introduce a new approach to model-free path-dependent 
option pricing. We first introduce a general duality result for linear opti¬ 
misation problems over signed measures introduced in Q and show how 
the the problem of model-free option pricing can be formulated in the new 
framework. We then introduce a model to solve the problem numerically 
when the only information provided is the market data of vanilla call or put 
option prices. Compared to the common approaches in the literature, e.g. 
14], the model does not require the marginal distributions of the stock price 
for different maturities. Though the experiments are carried out for simple 
path-dependent options on a single stock, the model is easy to generalise for 
multi-asset framework. 


1 Introduction 

In this paper we introduce a new approach on model-free pricing of discrete 
path-dependent options. There is a considerable amount of research done 
on the theory of the problem, but, to our knowledge, not much results are 
known for the numerical solution for the general case. This work is based 
on a general duality relation introduced in Q. This duality relation serves 
us twice in this work. First we show how the optimisation problem of find¬ 
ing model-free lower and upper bounds on pafh dependenf opfions can be 
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formulated in their framework, which allows one to get the dual problem 
straightaway from the primal. 

Next, we apply the so-called ’Constraints via Bounds on Integrals’ model to 
solve the optimisation problem numerically. This model allows one to con¬ 
struct constraints for marginal distributions, which can be obtained by cal¬ 
ibrating the empirical risk-neutral distributions. But instead we use market 
data of European call or put prices as constraints. This allows us to compute 
the tightest possible bounds in the model-free framework, i.e., given only the 
market information, any narrower bound has to consider constraints other 
than those provided by the market data and that the stock price process is a 
martingale. 


2 Problem Formulation and Duality 


Consider an exotic path-dependent option that depends on the value of a 
stock price at some discrete times ti <t 2 < ■ ■ ■ < t„. Let X = (Xi,... ,X„) be 
a n-dimensional random variable denoting stock prices at times (fi,... ,f„) 
respectively, and let h{Xi,. .. ,X„) be the payoff of the given option. 

In practice, in a no-arbitrage framework one assumes that the stock price 
follows some specific model, one calibrates this model to market data and 
finds the fair price of an option, i.e., find a risk neutral probability measure 
Q on M” and calculate the price as 



( 1 ) 


Under Q the stock price process „ is assumed to be a discrete mar¬ 

tingale. 

The problem of model-free option pricing is to find bounds on the prices of 
the option without assuming any particular model for the stock price process, 
i.e., find maximum and minimum values over all models that are consistent 
with the observed market data. Assuming knowledge of market call option 
prices, one can recover a distribution function of X, for every i = 

(see, e.g., [51), which will represent the marginal distributions of a risk neu¬ 
tral equivalent probability measure of X. In this framework, denoting by 
■ ■ )^«) the set of martingale measures with marginals 
the problem becomes 



( 2 ) 


We call the above problem the primal problem. 

For any i= define the set of Fi integrable functions from M ^ M by 

Mi and for any y = 1,..., n — 1 the set of bounded measurable functions from 
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—>■ M by S^j. For any hi G i= and gj G j = I,... ,n — I 

consider the function {xi,... ,x„) : M” —)• M defined by 


n n—l 

^{hl),{gj)ixi,. . . ,Xn) = Yj^i{Xi) + £ gj{xu. . . ,Xj){Xj+l -Xj). 
i=l j=l 


Consider the following optimisation problem: 

n 

sup £ 

S.t. <h{x\,...,Xn). 

We have the following theorem. 



( 3 ) 


( 4 ) 

( 5 ) 


Theorem 1 (Theorem 1, |T|) Assume ,..., are Borel probability mea¬ 
sures on R so that ,... ,^n) non-empty. Let /j: R” —>■ (—oo^oo] be 

a lower semi-continuous function so that 

h{xi,...,Xn) > -/:• (1 + |.ri| H-h|x„|), 

on R”/or some constant K. Then the duality gap between the primal problem 
(Ell and the dual (01 equals zero. Moreover, the primal optimal value is 
obtained, i.e. there exists a martingale measure ^ G ^ ,..., ^n) such 

that the optimal value of dUl is equal to 

A proof of the theorem can be found in [T] . 

In order to solve the dual problem numerically, we need to make some dis¬ 
cretisation of this infinite dimensional problem. One way of doing so (as 
is described in |T|) is to decompose the functions (/j,)” and (g;)” ^ over a 
finite dimensional basis. Since (g;)” ^ are assumed to be continuous, we 
can choose a polynomial basis, i.e. 

gj{xi ,... ,Xj) = Y^g)e){xi ,... ,Xj), 
k 

where e^j are functions of the basis, and g^j are the coordinates of gj with 
respect to the chosen basis. 

For functions (/z,)” we can choose call option payoffs as basis 

=Y^h'l{xi-k^)^. 

k 

With this approximation the dual problem reduces to a linear programming 
problem that can be solved via interior-point methods, the simplex algorithm 
or another algorithm of choice. 

In this work we take a different approach to the numerical solution of the 
problem of arbitrage-free pricing. The approach is based on the model de¬ 
scribed in Q, which shall be presented in the next section. 
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3 A General Duality Relation 


In this section we will see that the duality relation between (|2]) and (|4l) in¬ 
troduced by |T] in a special case of a general duality relation for mass trans¬ 
portation problems presented in Q, thus yielding an alternative derivation 
of®. 

First we will describe the duality results, and then will show how to obtain 
the duality results presented in the previous section using this framework. In 
this section we also present a model that lends itself to the numerical solu¬ 
tion of the resulting infinite-dimensional optimisation problems, which have 
special structure. 

Let (r, ©) and (L, S) be complete measure spaces, and let A : F x 

<I>—>-M, a: F^M, B: rx<I>^]R, f?: and h\ <I> —> M be bounded 

measurable functions on these spaces and the corresponding product spaces. 
Let and be the set of signed measures with finite variation on 

(<I>,5^), (F, ©) and (Z,S) respectively. 

Consider the following pair of optimisation problems over and x 
respectively. 



s.t. [ A{y,x)d^{x) < a{y), 

J<t> 

[ B{z,x)d^{x) = b{z), 

J<t> 

^> 0 , 

and 

[ a{y)dW{y)+ f b{z)dJy^{z), 

Jr Jr 

s.t. J^A{y,x)d'^{y) + J^B{z,x)dy{z) >h{x), (vG<I>), 

^>0. 


(V) 


The infinite-programming problems ® and ® are duals of each other. 

Theorem 2 (Weak Duality, Theorem 1, [O) For every ^-feasible measure 
and every ^-feasible pair we have 


/o 


[ h{x)d,^{x) < [ n(y)df#(y)-I- [ b{z)dy’{z). 

Jo Jr Jr 

Proof. Using Fubini’s Theorem, we have 

h{x)dJ^{x)< [ A{y,x)d{^ X ^){y,x) + [ B{z,x)d{y x .^){z, 

JrxO JExO 

< J^a{y)d'i^{y)+ J b{z)dy{z). 
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□ 

In various special cases, such as those discussed in EJi], strong duality is 
known to hold subject to regularity assumptions, that is, the optimal values 
of ® and ([T]) coincide. Another special case under which strong duality 
applies is when the measures ^ and 5^ have densities in appropriate 
Hilbert spaces, see the forthcoming DPhil thesis of the second author |'8]. 
Note also that the quantifiers in the constraints can be weakened if the set of 
allowable measures is restricted. For example, if is restricted to lie in a set 
of measures that are absolutely continuous with respect to a fixed measure 
% € then the quantifier (y € F) can be weakened to (%-almost all 

r). 


3.1 Application of the General Duality Result to Model 
Free Option Pricing 

In this subsection we show that the duality result for no-arbitrage bounds of 
path-dependent options can also be obtained in this framework. 

In the original Problem ([21) we have two types of constraints. The first type 
is the requirement of the joint probability measure ^ to have marginal dis¬ 
tributions =^ 1 ,... The second type of the constraint is for to be a 
martingale measure. 

The marginal constraints can be formulated as follows 

s-t- [ l{jc,<z}iz^x)d^{x) =Fi{z), (zGM,/= l,...,n) (8) 

Jm 

For the martingale measure constraints, observe, that is a martingale mea¬ 
sure if and only if for every k = 1,..., n — 1 and every Borel measurable set 


/ l^{xi,...,Xk)iXk+i -Xk)d^{xi,...,Xn) =0, 


(9) 


which is equivalent to (lT|, Lemma2.3): for every z,y € such that z<y 

/ ^z<xW<y{xi,---,Xk){xk+l-Xk)d^{xi,...,Xn)=0, ( 10 ) 

JR" “ “ 


where by x\k we denote the first k coordinates of x and the comparison of the 
vectors is component-wise. 

Now we are ready to formulate the Problem ^ in the new framework. 
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We aim to find the dual problem of the following optimisation problem. 


inf / 0(x)d^(x) 

JR" 

s.t. [ =Fi{t), (tGM,/GN„) 

JR" 

J^jz\k<xk<y\ki^^y^^)i^k+i -Xk)d^{x) = 0 , (y,zeW‘,k<n) 

^ > 0 . 

Denote <I) = M” and r = {l,...,2n — l}xMxM"xM". 

We eonstruet the funetion b{k,t,z,y) = Fk{t) ifk<n, and b{k,t,z,y) = 0 for 
k> n. 

Further, for k <n we have B{k,t ,z,y,x) = Ix^xo For fc=l,...,n — Iwe have 
B{k + n,t,z,y,x) = fz\k<x\,<yk{xk+i-Xk) 

Now we eonstruet the dual problem. 

As a objeetive funetion we have 

[b{k,t,z,y)dy{k,t,z,y) = [ Fk{t)d^k{t,z,y). 

Jl. “lJRxR"xR" 

For k = define the measures y'k on Borel sigma-algebra sueh 

that y'kiD) = F^k{F> X M” X M”) for every D G =^(M) . As for the marginal 
problem, the signed measures S^k being of finite variation, the funetions 
Skit) = ^ii—°°it]) and the limits Sk = limf^oo5'<:(t) = =5^((—oo,-|-oo)) are 
well defined and finite. Furthermore, using Imrt^-ooFkit) = 0 and 
lim(^+ooiy;(t) = 1, we have 

£ f Fkit)dS^kit) = £ (f,(05, (Oil: - / 5,(t)dF,(o) 

= / Skit)dFkit) 

= L / isk-Skit))dFkit), 

and likewise for eonstraints. Writing /i,(t) = Sk — Skit), (k = the 

objeetive funetion beeomes 


£ / hkit)dFkit), 

ITTTl 
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and the left hand side of the constraints become 


B{k,t,z,y,x)Ay{k,t,z,y) 


^z\kMk<y\k^^k+\ -Xk)dyn+k{t,z,y) 
n—i 

= Y^hk{xk) + Y^ gk{xi,... ,Xk){xk+\ -Xk), 
k=i 


k=\ 
n—l f 

v/rxR"x 


t:=l 


where ... ,Xk) = /r«xR" K\k<x\k<y\k^^>^+^ -Xk)Ayk+n{z,y), k = l^...,n- 
1. Measures y'n+k, k = I,... ,n—l are defined such that ^k+n(D) = ^k+n (1^ x 
D) for every D G ,^(M” x M"). 

We obviously have that the functions g^, k = 1,... ,n — 1 are bounded. So 
we get the following dual formulation of our primal problem. 


sup y / hk{T:)dFk{T) 


( 11 ) 


n—\ 


s.t. Y,^kixk) + Y,Sk{xi,...,Xk){xk+i-Xk) <h{x), (.tGM”). 

k=0 k=\ 

Since the decision variables gi, ... ,gn-i have zero coefficient in the objec¬ 
tive function, we get that the dual problem (fTTI) is equivalent to 

sup \ Ef, [hk] :3gu- --gn-i s.t. < 0 L (12) 


. k=\ 


where 


n—l 


= Y^^kixk) + Y, Skixi,. ■ . ,Xk){Xk+l -Xk) 
k=r) k=\ 

and EFt{hk{t)) = j^h{t)dFk{t). 

Note that ^Y2\ is exactly the same dual as (|4l). 


3.2 Constraints via Bounds on Integrals 

In this section we present a particular case of the general duality result dis¬ 
cussed at the beginning of this section, specifically, where we have finifely 
many consfrainfs. In fhe nexf section we show how fo find no-arbifrage 
bounds of fhe opfions wifh piece-wise linear payoffs using fhis model. 
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Let us assume that <I> is decomposed into a partition <I> = Ijf^j E, of polyhe- 
dra E; with nonempty interior. Each polyhedron has a primal description in 
terms of generators, 

E; = conv q'n-) + cone (r'l,..., ) 

where conv(( 7 j,... ,q'^.) is the poly tope with vertices € M”, and 

cone(ri,...,r^.) = < £ UC : > 0 (m G N^,) 

is the polyhedral cone with recession directions G M”. Each polyhedron 
also has a dual description in terms of linear inequalities, 

k 

E,- = n e M" : {fi,x) > 1)] , 

;=i 


for some vectors f \ G M” and bounds £'■ G M. The main case of interest is 
where E,- is either a finite or infinite box in M” with faces parallel to the 
coordinate axes, or an intersection of such a box with a linear half space, in 
which case it is easy to pass between the primal and dual descriptions. Note 
however that the dual description is preferable, as the description of a box in 
M” requires only 2n linear inequalities, while the primal description requires 
2" extreme vertices. 

Consider the problem 

sup f h{x)d£^{x) (13) 

•f't’ 

s.t. [ (l)s{x)d^{x) <as, {s = l,...,M), 
s.t. [ \j/t{x)d^{x) =bt, (f = l,...,N), 

JO 

[ 1 d ^ (x) = 1, 

v/o 

^>0, 


where the test functions \j/t ^c piecewise linear on the partition <I> = E/, 

and where —h{x) and the test functions are piecewise linear on the infinite 
polyhedra of the partition, and either jointly linear, concave, or convex on 
the finite polyhedra (i.e., polytopes) of the partition. The dual of (fTSl) is 


inf 

(>',z)eK'^+"+> 


M N 


s=\ t=\ 


( 14 ) 


M N 

S.t. ^y^^,(x) + ^z,v^i(x)+zoloW-/iW >0, (x€<I>), 

i=i t=\ 


y >0. 
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Note that (fT^ is a semi-infinite programming problem with infinitely 
many variables and finitely many eonstraints, while (fTdl) is a semi-infinite 
programming problem with finitely many variables and infinitely many con¬ 
straints. However, the constraint of the problem (fTdl) can be rewritten as 
co-positivity requirements over the polyhedra E,, 

M N 

Y^ys(j)s{x) + Y,ZtVt{x)+Zo'^^{x)-h{x)>0, (vGS,-), {i=l,...,k). 

5'—1 t—l 

Next we will see how these co-positivity constraints can be handled numer¬ 
ically, often by relaxing all but finitely many constraints. 

In what follows, we will use the notation 

M N 

%,z{x) = + Y^ZtWt{x) +Z0 - h{x)- 

i=i t=\ 


3.2.1 Piecewise Linear Test Functions 

The first case we discuss is when and h\zj are jointly linear. Since we 
furthermore assumed that the functions Yt\zi are linear, there exist vectors 
vj € M”, wj G M", g' G M” and constants G M, d} G M and e' G M such that 

(l>As,{x) = (v',v)+c', 

Vtkiix) = {w\,x)+dl, 
h\zi{x) = {g\x) + e\ 


The copositivity condition 

M N 

Y,ys(l>s{x) + Y,^tVtix) + Zo'^<S>{x) - h{x) >0, (xGS;) 

s=l /=! 

can then be written as 


{fvx)>i\ (; = 1,...,A:0 


M 


N 


M 


L ysvi + £ Ztw] - g' , X ) > e' - £ y,c' - £ Ztdi - zo 


\.s=l / s=l 

By Farkas’ Lemma, this is equivalent to the constraints 


N 

L 

t=i 


M N ki 

£ ysv's + £ ztw\ - g' = £ X‘jfj, 

s=l t=\ 7=1 

M N ki 

«^' - £ - £ 1,4 - ZO < £ 

i=l f=l 7=1 

Aj> 0 , U = h---,ki), 


(15) 

(16) 
(17) 
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where Aj are additional auxiliary decision variables. 

Thus, if all test functions are linear on all polyhedral pieces S/, then the dual 
(D) can be solved as a linear programming problem with M+N+\ h 

variables and k{n + 1) linear constraints, plus bound constraints on y and the 
Xj. More generally, if some but not all polyhedra correspond to jointly lin¬ 
ear test function pieces, then jointly linear pieces can be treated as discussed 
above, while other pieces can be treated as discussed below. Details of effec¬ 
tive solution methods of the above constructed linear optimisation problem 
can be found in l3i|. 

In the next section we show that the problem of model-free option pricing 
with piece-wise linear payoffs can be approximated in the framework of this 
model. The numerical experiments we present in Section 0 show that in 
some cases the model solves the problem exactly. 

4 Numerical Implementation 

In practice we are not given European option prices for every strike, hence 
we cannot have exact knowledge of marginal distributions. Moreover, for 
some maturities the number of traded liquid products can be so small that 
any interpolation of option prices will contain additional errors. To avoid 
introducing additional estimation errors in the model, and to use only the in¬ 
formation given by the market, we replace constraints for the marginal distri¬ 
butions by constraints which represent calculation of vanilla option values. 
Taking into consideration the above observations, we construct the following 
problem 

sup [ h{x)d^{x) 

^ JR" 


s.t. / {Xi-kl)^d^{s) 
JR" 

<c/, / = i.. 

.,n, j = l,.. 

■^mj 

f {xi — kj)^d^{s) 

>C/, / = !.. 

.,n, j=l,.. 

■^mj 


[ (l)i{s)d^{s) <di, i=l,...,N^ 

Jr" 

[ (l)i{s)d^{s)>5i, i = l,...,N^ 

jR" 

[ ^z<s\k<y{z,y,^){xk+i-Xk)d^{s) =0, (y,z G < n) (18) 
JR" “ “ 

[ \d^{s) = \, 
jR" 

^> 0 , 

where are the strikes of available in the market European call op¬ 
tions, {CjyjLi and are their respective market bounds of prices (bid- 
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ask prices) and m, is the number of such options. 

Except of call or put options we might observe other options in the mar¬ 
ket. Let them have payoff functions ^i{s) and bid-ask prices 5, and Si, 
i = For each of these additional observed options we added 2 

more inequality constraints. 

Note that the constraints representing call option prices fit our model (the 
functions are piecewise linear). To render the martingale constraints com¬ 
patible with our model, we approximate them by requiring that (fTOl) hold for 
only finitely many boxes in M” of our choosing. 

We impose a lower and upper bound on the stock price X, for each i = 
0 will represent a natural lower limit, and for the upper limit, we 
can take a large enough value that the stock price is statistically unlikely to 
exceed. We also make sure the strikes are included in the set of points in our 
discretisation. 

Let 0 = d}<d}<---< df‘ be the set of discretisation points for X,-. Denote 
<I> := {5 : d/ < X/ < rf”', (/ = 1,...,n)}. 

Then, the martingale constraint for k = 1 is 



1 


{xeR" <Xi <d/,{i=\,...,k)} 


[s){xk+i -Xk)d,^{s)=0, 


i^<ji<^i<ni,\/i). 


For each k we thus have Mk := IlLi “1) equality constraints, which are 
all represented as integral inequalities with piecewise linear functions. 

We saw that after discretisation all our constraints fit the model described in 
Section [T2] Further, if the payoff functions h{s) and ^i{s), / = !,... ,N,j, 

are piecewise linear (put/call or binary barrier options, Asian options, look- 
back options etc.), then we can solve the dual problem as described in Sec¬ 
tion |T2j If the payoff function h{s) is concave (as for options on realised 
variance), then we can find the upper bound with our method, and if it is 
convex, we find the lower bound. For all other cases we have to make ap¬ 
proximations, for example, first order Taylor approximation. 

Now we can write the dual problem as follows: 


inf 

y,y,z,z,t,w 


S.t. 




L L (yi cj -yJi^i) + Y. ^ 

i=lj=l (=1 

n Hi ^(p 

E Ee; -yl)(x.-Kiy E(s-s.)ft(»)+ 

/—I 7—1 i —1 

n— 1 

L L K.Mk<Pr+,i^)i^k+l-Xk)+W>h{s), 
k=l tg4 


Vi' G <I>. 


Interested readers can find efficient algorithms to solve the above de¬ 
scribed optimisation problems in the DPhil thesis of the second author [8]. 
Since the size of the linear optimisation problems grows rapidly with the 
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dimension, memory cost becomes an important issue. For this situation it is 
useful to implement the delayed column generation method of the simplex 
algorithm. Instead of creating the whole matrix of the problem, the delayed 
column generation algorithm gradually adds columns to the tableau while 
the simplex algorithm takes iterates. In practice, one does not need to con¬ 
sider all the columns in order to achieve the optimal solution to a desired 
accuracy. Moreover, since a large number of iterations works on just a small 
number of columns, the method gains a great advantage in saving time by 
not using the secondary memory. 

Arbitrage-free pricing is often used as a means to identify arbitrage oppor¬ 
tunities. Suppose some path-dependent option with payoff h{s) is traded in 
the market with price h!^. By implementing our model for this option we 
find no-arbitrage bounds [L,U]. If ^ [L,U], then there is an arbitrage 
in the market (the dual problem also provides with the strategy to make 
risk-free money in this case, but we will not go into details in this work, 
interested reader can refer to idl). Hence, often we are interested in finding 
ouf whefher, for a given h^, L < <U. For fhis particular purpose, fhe 
defailed descripfion of fhe implemenfafion of Nesferov’s algorifhm fo linear 
opfimisafion problem can be found in fS ]. 


Remark 3 Note that one can introduce additional constraints into the model 
depending on additional assumptions one wishes to make on the behaviour 
of the stock. For example, it is easy to include lower and/or upper bound on 
the volatility of the stock. 


5 Payoffs as Piece-wise Linear Functions 

In fhis secfion we describe differenl payoffs fhaf are piece-wise linear func- 
fions. Obviously, call and pul oplions hi lo fhis calegory. 

5.1 Barrier Options 

For given lower and upper barriers B\ < B 2 (if the option is defined for one 
barrier only, e.g. down and ouf, up and in, Ihen we suppose fhe olher barrier 
is —00 or 00 depending on fhe conlexl), barrier opfions have payoff 

h{s)=H{X)fl\B,<x,<B„ 

1=1 

where H{X) is fhe final payoff. For digifal barrier opfions we have H{S) = 
1, 'is and for call and pul barrier oplions wilh slrike K we have 
H{S) = (A„ —and//(S') = {K — X„)^ respectively. 

If is obvious fhaf fhese opfions have piece-wise linear payoffs. 
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5.2 Lookback and Asian Options 

Lookback options have payoffs that depend on the X^ax = maxjXi,... ,X„} 
and/or X^in = min{Xi ,...,X„}. Asian options are defined as having payoff 
functions at maturity which depend on the average value Xa = of 

the stock price over the path. Table [T] describes payoffs for different look- 
back and Asian options. The following lemma proves that all payoffs de- 

Table 1: Lookback and Asian Option Payoffs 



Call 

Put 

Lookback Fixed 

msix{X^ax-K,0) 

max{K - X,nin,0) 

Lookback Floating 

^min 

Xmax Xfi 

Asian Fixed 

max(2C4 — K, 0) 

max{K — XA,0) 

Asian Floating 

Xu-Xa 

Xa-X„ 


scribed in the Table [T] are piece-wise linear functions on <I>. 

Lemma 4 Let ^ be a polytope in M” and be further divided into finitely 
many polytopes <I> = (U^i nO- Let function / : <I> —^ M and g : <I> —>■ M /jc 
piece-wise linear on these polytopes. Then, the functions max{/(x),g(v)} 
and min{/(v:),g(x)} are also piece-wise linear on some finite partitioning 
of^ = {\jt,nt). 

6 Numerical Experiments 

In this section we present numerical results for barrier options. We con¬ 
sider 2 and 3 dimensional cases. Numerical solutions of the linear program¬ 
ming problem are carried out with interior point methods using SeDuMi 
and SDPT3 solvers on MATLAB. For higher-dimensional or finer discreti¬ 
sation problems it is better to use Nesterov’s smoothening technique for non¬ 
smooth functions, IZlii. 

First we obtain our numerical results for simulated market data, i.e. we 
compute call option prices with Black-Scholes formulae for different time 
horizons. Numerical experiments on simulated market data is of particular 
interest because in this case we can also compute the true prices of our ob¬ 
jective derivative (i.e. the one that we aim to find upper and lower bounds 
on). Hence, we can compare the price of the derivative in Black-Scholes 
model with the lower and upper bounds obtained by our model. 

Tables |2] and [3] show numerical results for different barrier options for 2 
and 3 time steps respectively. To obtain these results we only used market 
given strikes as discretisation points, i.e. for all / = l,...,n and j = I,... ,ni. 
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Lower 

Upper 

Black-Scholes 

Double no-touch digital 

0.282622 

0.612447 

0.4232 

Double no-touch call 

0 

0.527483 

0.202 

Variance Swap 

0 

0.187674 

0.0919 


Table 2: Lower and Upper Bounds on 2-time step barrier options compared to 
Black Scholes value. 



Lower 

Upper 

Black-Scholes 

Double no-touch digital 

0.0610184 

0.533453 

0.2732 

Double no-touch call 

0 

0.43506 

0.1 

Variance Swap 

0 

0.326711 

0.1369 


Table 3: Lower and Upper Bounds on 3-time step barrier options compared to 
Black Scholes value. 



Barrier Call 

Barrier Digital 

Variance Swap 

2 

0.554805 

0.90002 

1.21195 

3 

0.485545 

0.855626 

1.24816 


Table 4: Lower bounds on variance swap, double no-touch barrier digital and cal 
options for 2 and 3 time-steps. The results are obtained using market data for 
Apple stock. The stock price is 115 and the barriers are 95 and 125. 
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dj = kj. 

We carried out the experiments for simulated option prices with spot price 
50, volatility 30% and time horizons with time steps 6 months. Strikes are 
chosen at all even numbers between 30 and 60 inclusive, while lower and 
upper bounds were chosen as 34 and 56 respectively. 

In Table|4]we present lower bounds on double no-touch call, digital and vari¬ 
ance swap options using market data for call options on Apple stock (APPL). 
The spot price is 115 and the barriers considered are 95 and 125. We used 
call option prices for maturities 20/12/14, 20/02/15 and 17/04/15. The 
calculation date is 20 November 2014. 

Note that the third option considered here is the variance swap option, which 
has neither piecewise linear nor piecewise concave payoff. In this case, we 
approximate the payoff function by a piecewise linear function with a linear 
restriction over each polytope. 

After carrying out some numerical experiments, we observe that as long as 
prices of call options for the strikes that coincide with barriers of our options 
are available, any finer discretisation than the one described above leads to 
no further tightening of the bounds, that is, the solution of the optimisation 
problem remains the same. This observation, a proof of which is the sub¬ 
ject of further research, implies that our method yields the tightest possible 
bounds that are compatible with the given market data. This is to say that for 
any other bound I > L and u < U one can construct a martingale measure 
for stock prices such that under this measure we get prices for our barrier 
options that are out of these new bounds I and u, where L and U are the 
lower and upper bounds obtained by our model. 

Taking into consideration the observation above, in the case when we don’t 
have barriers as strikes in our market data, we should make some interpola¬ 
tion at the beginning to find approximafe call option prices for fhe barriers 
as sfrikes and fhen consfrucf our model. This approach is heller lhan solving 
fhe problem by making a finer approximation. 


Remark 5 Note that in practice information on call option prices provided 
by the market and the martingale constraint produce very wide bounds, so 
practically one will never observe an arbitrage, i.e. the price of the objective 
option will always be within the bounds. But as can be observed in the 
results that use call option prices obtained from Black-Scholes model, the 
real value of the option is approximately the median of the lower and upper 
bounds. 
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7 Conclusions 


In this work we presented how the general duality results from ||3l can be 
applied to model-free path-dependent option pricing. We also introduced 
the so-called ’constraints via integrals of piece-wise concave functions’ to 
solve the resulting semi-infinite dimensional optimisation problems numeri¬ 
cally. Moreover, our numerical experiments provide strong evidence that the 
bounds obtained via our method are the tightest possible bounds compatible 
with the observed market prices and the no-arbitrage assumption. 
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